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Evidence for the same universal behavior of 2d Yukawa (Y2) and Gross-Neveu (GN) models in a certain range 
of couplings, particularly also for k < 0, is presented. 
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1. Introduction 

Yukawa models in 2d (Y2) with chiral Z(2) and 
U(l) symmetry are investigated on the lattice. 
They are natural extensions of the GN models, 
which are specially interesting because of their 
asymptotic freedom and the dynamical mass ge- 
neration. Some years ago the equivalence bet- 
ween Y2 and GN models has been suggested in 
[1] using a mean field (MF) analysis. Recently 
we have found numerical evidence for the asym- 
ptotic freedom of the Y2 models in a large range 
of couplings [2,3]. This property is suggested in 
the effective potential approach [4] and by our 
MF method [3], which is reliable in 2d due to 
infrared divergence of the fermion loop inducing 
long-range effective scalar field interaction. 

In this paper we discuss the question whether 
Y2 and GN models are in the same universality 
class for a certain range of couplings. Most of 
the results were obtained by simulating the action 
with staggered fermions and hypercubic Yukawa 
coupling, e.g. in the Z(2) case: 
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where k is the hopping parameter, A the quartic 
selfcoupling of the scalar field, ry^ X the sign factor 
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Figure 1. Sketch of the phase diagram of Y2 mo- 
dels with hypercubic Yukawa coupling. 

for staggered fermions and <f)^ the average of <f) on 
the hypercube with the lower left corner x. The 
number of fermions in the continuum is Np = 
AN, a = I..N labeling the fermion species. 

In this representation the GN model arises in 
the limit k = A = 0, while at A = cxd, t/ = eq. 
(1) describes the Ising or the XY model, Kc(X) 
being the critical line. 

2. Properties of Y2 models 

Works on the GN models and on the Y2 models 
[1-4] lead us to the phase diagramm sketched in 
fig. 1. It is based on l/N expansions, which can 
be justified only for A = 0(l/N), MF approxima- 
tions, which can give very good results in 2d [5], 
and numerical simulations of our group. 

Due to the logarithmic divergence in 2d, the 
phase diagramm differs from that of higher di- 
mensions. At small y there is no paramagne- 
tic phase in the Z(2) case, the whole sheet 
{t/ = 0,K < Kc(X)} being critical. At k —Kc(X) 
another critical surface separates the ferromagne- 
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tic (FM) from the ferrimagnetic (FI) phase. In 
the U(l) case the phase structure is similar, we 
have above —Kc(X) the spin wave (SW), below the 
staggered spin wave phase. 

Fermion mass is generated dynamically in both 
Y2 models [3,5,6] for arbitrarily small y. The con- 
tinuum limit can be performed by approaching 
the critical surface at j/ = 0, k < Kc(X). 

Y2 models are asymptotically free like the GN 
models [7,2,3]. At fixed k and A, for k < Kc(X) 
the MF approximation [3,5] suggests in the con- 
tinuum limit the scaling law: 

amp - exp ( -77^ , \\ 2 ) ' (2) 

x(k, A) being the susceptibility in the scalar theo- 
ries (y = 0) and /3o the first coefficient of the 
/3-function in the GN models. For k > it is 
X = Z^/a'^m'^, but the MF prediction holds also 
for K < 0. 

3. Universality 

The phase structure in fig. 1 suggests a classi- 
fication of the possible universality classes in the 
Y2 models, depending on where and how the con- 
tinuum limit is performed. At j/ = we have the 
usual theories, the phase transition at Kc(\) 
is for A > Ising like in the Z(2), Kosterlitz- 
Thouless (XY) like in the U(l) case. The Gaus- 
sian fixed point is at A = 0, k = 1/4. 

When approaching the critical line Kc(X) from 
positive y, both, fermions and bosons will remain 
in the spectrum and the situation is very inte- 
resting. The mass ratio of fermions and scalar 
particles must depend on the way the continuum 
limit is performed, characterizing a whole family 
of universality classes, similar to the case descri- 
bed by A. Sokal at this conference [8]. 

At the critical surface around k ~ —Kc(X), 
y > 0, the fermion mass doesn't scale. In the 
corresponding continuum limit only the stagge- 
red scalars would remain in the spectrum. 

In the continuum limit performed on the cri- 
tical surface y = 0, k < Kc(X) only the fermions 
and their bound states survive in the spectrum. 
Fermions scale like in the GN model, but it is 
not clear whether GN and Y2 models belong to 



the same universality class in this region. Two 
aspects of this question will be discussed in the 
following: whether /3o of Y2 and GN models is 
the same, and whether the mass ratio of the first 
two states is the same. 

The first coefficient of the /3 function can be 
interpreted as a critical exponent. To compare 
it with that of the GN model one has to consi- 
der the GN scaling law: amp ~ exp ^"2^3^)' 
with g being the usual 4 fermion coupling. The 
l3o coefficients would be the same if the effective 
scaling variable in the Y2 case is g = ^JxV- This 
is indeed suggested by perturbation theory in y. 
The fermion 4-point function is modified already 
on the tree level in y by the scalar field. At zero 
momentum it is: 

= 0) = j/2G(f)(K = 0) = y\{K, A). (3) 

As eq. (2) is valid only in infinite volume at 
very small j/, we have to modify it in order to 
compare it with finite lattice simulations [3]. We 
use lattice sums instead of momentum integrals 
and the full response function /(-ff) determined 
in numerical simulations in order to make pre- 
dictions about amp with the MF method. The 
selfconsistency eq. 

= /[#j/amf ^(^sm2p^ + a2^^)-i]^ (4) 

with amp = y{<^'), must be solved numerically, by 
recursion. The result is compared to amp . But as 
this method is only an approximation compara- 
ble to 1 fermion loop results, we have to consider 
eventual corrections to it. Therefore we intro- 
duce a parameter which quantifies the deviation 
from the MF prediction by taking amp = r • t/((/)) 
and fitting several amp values with the function 
obtained by solving (4). Fig. 2 shows a typical 
fit and the resulting values of r in such fits are 
compiled in table 1. 

It is clear that r k, \ means that the scaling 
oi amp is GN like. But certain deviations can be 
expected, because MF is not perfect. An estimate 
for these deviations is given by the fit in the GN 
case, where r differs from 1 by about 5%. Such 
a deviation for other k, A values thus does not 
contradict the hypothesis, that /3o is the same as 
in the GN case. 
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Figure 2. Scaling behavior of amp for k = —0.2 
and A = cxD on 16^, 32^ and 64^ lattices. The lines 
are fits for the parameter r defined in the text. 



We decided to calculate the r values for A = cxd, 
at maximal bare scalar selfcoupling. The results 
for K = —0.4 and k = —0.2 show that r 1 
whithin the limits allowed by the GN value at k = 
A = 0. This strongly supports the assumption 
that Y2 and GN models are equivalent even at 
negative k. 



K A 


162 322 




-0.4 CXD 
-0.2 CXD 
CXD 
0.2 CXD 


0.957(3) 0.956(5) 
0.971(5) 0.970(7) 
0.955(5) 0.95(1) 0.94(3) 
0.94(1) 0.93(2) 0.93(4) 
0.91(4) 0.86(4) 0.85(6) 



Table 1: The values of the parameter r. 

At K = 0.2 the parameter r differs considera- 
bly from 1. The quality of the fits was poor, 
i.e. the data could not be described well by the 
MF approach. This is probably an effect of the 
crossover to the universality classes at k = Kc(X). 
It could also explain the relatively large deviation 
(r ^ 0.93) at K = 0, A = cxD. 

The mass ratios are universal quantities and 
should be the same if both models are in the same 
universality class. We tried to measure the first 
one. In fig. 3 m-pp/niF is plotted against another 
dimensionless quantity, mpL. Data from all lat- 
tices with L/ a = 8, 16, 32, 64 is plotted with the 
same symbols, so the mass ratios seem to scat- 
ter more than they in fact do on a given lattice. 
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Figure 3. The mass ratio m-pp/niF vs. mpL. 

All K values give the same results and they are 
consistent with those of the GN model (circles). 

Concluding we can say that we found substan- 
tial evidence for the assumption that Y2 and 
GN models are in the same universality class 
for K < Kc(X). Particularly interesting is the 
fact that this holds also for negative k, which is 
beyond the scope of the considerations in conti- 
nuum [1]. 

In 2d Yukawa theories reduce to the simpler 
GN models. The situation is different from that 
in 4d, where the Yukawa theories are those with 
the right number of parameters, not the 4 fermion 
(Nambu - Jona-Lasinio) models [9]. 
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